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The properties of jets produced inp + p collisions at y's=200 GcV arc measured using the method 
of two particle correlations. The trigger particle is a leading particle from a large transverse momen- 
tum jet while the associated particle comes from either the same jet or the away-side jet. Analysis of 
the angular width of the near-side peak in the correlation function determines the jet fragmentation 
transverse momentum jT- The extracted value, = 585± 6(stat) ± 15(sys) MoV/c, is constant 

with respect to the trigger particle transverse momentum, and comparable to the previous lower ^/s 
measurements. The width of the away-side peak is shown to be a convolution of jr with the frag- 
mentation variable, z, and the partonic transverse momentum, /ct. The {z) is determined through 
a combined analysis of the measured tt" inclusive and associated spectra using jet fragmentation 
functions measured in e^e~ collisions. The final extracted values of kr are then determined to also 
be independent of the trigger particle transverse momentum, over the range measured, with value 



of = 2.68 ± 0.07(stat) ± 0.34(sys) GeV/c. 
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I. INTRODUCTION 

The goal of this paper is to explore the systematics of 
jet production and fragmentation in p + p collisions at 
y/s=200 GeV by the method of two-particle azimuthal 
correlations. Knowledge of the jet-fragmentation pro- 
cess is useful not only as a reference measurement for a 
similar analysis in Au + Au collisions, but can be used 
as a stringent test of perturbativc QCD (pQCD) calcu- 
lations beyond leading order. 

The two-particle azimuthal correlations method 
worked well at ISR energies (■y/s=63 GeV) and below 
0,1^1^, where it is difficult to directly reconstruct jets, 
but has not been attempted at higher values of y/s. 
This method is also suitable for jet-analysis in heavy 
ion data where the large particle multiplicity severely 
interferes with direct jet reconstruction. 

With the beginning of RHIC operation, heavy-ion 
physics entered a new regime, where pQCD phenomena 
can be fully explored. High-energy partons materializ- 
ing into hadronic jets can be used as sensitive probes of 
the early stage of heavy ion collisions. Measurements 
carried out during the first three years of RHIC oper- 
ation at ■y/SArAr=130 and 200 GeV exhibit many new 
and interesting features. The high-p^ particle yield was 
found to be strongly suppressed in Au -\- Au central 
collisions 0. Furthermore, the non-suppression of the 
high-pT particle yield in d -I- Au induced collisions @ 
confirmed that the suppression can be fully attributed 
to the final state interaction of high-energy partons with 
an extremely opaque nuclear medium formed in Au+Au 
collisions at RHIC. 

Other striking features found in RHIC data are the 
large asymmetry of particle azimuthal distributions 
which is attributed to sizable elliptic flow Q and the 
observation of the apparent disappearance of the back- 
to-back jet correlation in central Au + Au collisions :S]. 

Many of the above mentioned observations can be 
explained by a large opacity of the medium produced 
in central Au + Au collisions which causes the scat- 
tered partons to lose energy via coherent (Landau- 
Pomeranchuk-Migdal 0) gluon brcmsstrahlung 10, 
[11) . It is expected that the medium effect will cause 
the apparent modification of fundamental properties 
of hard-scattering like broadening of intrinsic parton 
transverse momentum kx |l2l IT^ and modification of 
jet fragmentation |l4j . Thus the measurement of jet 
fragmentation properties and intrinsic parton trans- 
verse momentum fc^ for P + P collisions presented here 
provides a baseline for comparison to the results in 
heavy ion collisions, helping to disentangle the complex 
processes of propagation and possible fragmentation of 
partons within the excited nuclear medium. 

This paper is organized as follows: Section ^ dis- 



cusses the method of two-particle correlations and the 
relations between jet properties and the angular cor- 
relation between parton fragments. The details of the 
PHENIX experiment relevant to this analysis are out- 
lined in section IhH Section Hvl deals with the analysis 
of the correlation functions extracted from the p -\- p 
data and an evaluation of the (jx) and (kx) quantities. 
The combined analysis of the inclusive and associated 
PT-distributions is discussed in section Ivl and the sen- 
sitivity of the associated pT-distributions to the frag- 
mentation function is discussed in section IVll Section 
IVIII presents the resulting values of the partonic trans- 
verse momenta kx corrected for the mean momentum 
fraction (zt). Section [VIIII summarizes the results from 
this paper. 



II. JET ANGULAR CORRELATIONS 

Jets are produced in the hard scattering of two par- 
tons [13, lia, UM ■ The overall p + p hard-scattering 
cross section in "leading logarithm" pQCD is the sum 
over parton reactions a + b ^ c + d (e.g. g + q ^ g + q) 
at parton-parton ccnter-of-mass (cm.) energy \/|, 



dxidx2dcoa 6* 



2xiX2 



* Deceased 
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(1) 

where fa{xi), fb(x2), are parton distribution functions, 
the differential probabilities for partons a and h to carry 
momentum fractions xi and X2 of their respective pro- 
tons (e.g. w(x2)), and where 9* is the scattering angle 
in the parton-parton cm. system. The parton-parton 
cm. energy squared is s = X1X2S, where \/s is the cm. 
energy of the p + p collision. The parton-parton cm. 
system moves with rapidity y = (1/2) ln(xi/a;2) in the 
p -\- p cm. system. 

Equation^gives the px spectrum of outgoing parton 
c (emitted at 6*), which then fragments into hadrons, 
e.g. a TT*^. The fragmentation function (2,/i^) is 
the probability for a 7r° to carry a fraction z = p'" /p'^ 
of the momentum of outgoing parton c. Equation ^ 
must be summed over all subprocesses leading to a tt" 
in the final state. The parameter is an unphysical 
"factorization" scale introduced to account for coUinear 
singularities in the structure and fragmentation func- 
tions 0, |23| , which will be ignored for the purposes of 
this paper. 

In this formulation, fa{xi), fb(x2) and D'^ {z) repre- 
sent the "long-distance phenomena" to be determined 
by experiment; while the characteristic subprocess an- 
gular distributions, S"''(cos 6**), and the couphng con- 
stant, as{Q'^) = 25 inf^o^/ A ^' fundamental predic- 
tions of QCD 122, 122, [23] for the short-distance, large- 
Q^, phenomena. The momentum scale ~ p^ for 
the scattering subprocess, while ~ s for a Comp- 
ton or annihilation subprocess, but the exact meaning 
of tends to be treated as a parameter rather than a 
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FIG. 1: (color online) (a) Schematic view of a hard scat- 
tering event in the plane perpendicular to the beam. Two 
scattered partons with transverse momenta pr in the par- 
tons' center of mass frame are seen in the laboratory frame 
to have a momenta pTt and Ptb,- The net pair transverse 
momentum p^pair corresponds to the sum of the fey-vectors 
of the two colliding partons. The trigger and associated 
jet fragments producing high-pr particles are labeled as pTt 
and pTa.- The projection of fey perpendicular to pTt is la- 
beled as kxy The transverse momentum component of the 
away-side particle pVa perpendicular to trigger particle prt 
is labeled as pout. (b) The same schematics as in (sl), but 
the jet fragmentation transverse momentum component jxy 
of the trigger jet is also shown. 



dynamical quantity. 

Figure ^ shows a schematic view of a hard-scattering 
event. The transverse momentum of the outgoing scat- 
tered parton is: 



PT ^Pt 



sm( 



(2) 



The two scattered partons propagate nearly back-to- 
back in azimuth from the collision point and frag- 
ment into the jet-like spray of final state particles (see 
Fig. n^a) where only one fragment of each parton is 
shown) . 

It was originally thought that parton collisions were 
collinear with the p + p collision axis so that the 
two emerging partons would have the same magnitude 
of transverse momenta pointing opposite in azimuth. 
However, it was found 3] that each of the partons 
carries initial transverse momentum kx, originally de- 
scribed as "intrinsic" This results in a momen- 
tum imbalance (the partons' px are not equal) and an 
acoplanarity (the transverse momentum of one jet does 
not lie in the plane determined by the transverse mo- 
mentum of the second jet and the beam axes). The jets 
are non- collinear having a net transverse momentum 

It is important to emphasize that the (kr) denotes the 
effective magnitude of the apparent transverse momen- 
tum of each colliding parton. The net transverse mo- 



mentum of the outgoing parton-pair is a/2 • (fcr)- The 
naive expectation for the pure intrinsic parton trans- 
verse momentum based on nucleon constituent quark 
mass is about w 300 MeV/c [H H^]. However, the 
measurement of net transverse momenta of diphotons, 
dileptons or dijets over a wide range of center-of-mass 
energies gives (kr) as large as 5 GeV/c H^. It is com- 
mon to think of the net transverse momentum of a dilep- 
ton or dijet pair as composed of 3 components: 



pair 



T/ 



{kl 



t)- ■ - 

/ intrinsic 



/ soft ' 



t/nlo ' 
. (3) 

where the intrinsic part refers to the possible "fermi mo- 
tion" of the confined quarks or gluons inside a proton, 
the NLO part refers to the power law tail at large values 
of PTpair due to the radiation of an initial state or final 
state hard gluon, which is divergent as the momentum 
of the radiated gluon goes to zero, and the soft part 
refers to the actual Gaussian-like distribution observed 
as pt ^ 0, which is explained by resummation |27| . 

In the discussion below we will assume that the 
two components of the vector /ct, ^Tx and kxy are 
Gaussian distributed with equal standard deviations 
ciparton.id, in which case k'^=k'^^+k'^y is distributed 
according to a 2-dimensional (2D) Gaussian 26]. For 
the net transverse momentum of the jet pair, {Pr) pax-c ~ 



2partons,2d 



2rr2 



Note that the principal dif- 
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ference between the 1 and 2 dimensional Gaussians is 
that (fcTx) = (fcTy)=0, while {kr) ^ since kr is a 2D 
radius vector. 

The two components of fcy result in different exper- 
imentally measurable effects, fcyy leads to the acopla- 
narity of the dijet pair while /ctx makes the momenta 
of the jets unequal which results in the smearing of the 
steeply falling px spectrum. This causes the measured 
inclusive jet or single particle cross section to be larger 
than the pQCD value given by Eq.^ This was observed 
in the original discovery of high pT particle production 
at the CERN ISR in 1972 jH and led to much confusion 
until the existence and effects of k^ were understood. 

Before the advent of QCD, the invariant cross sec- 
tion for the hard-scattering of the electrically charged 
partons of deeply inelastic scattering was predicted for 
p + p collisions to follow a general scaling form: |i29i. i30| 



d^a 1 , , 
a P Pt VS 



(4) 



where xt = ^PtI \fs- The cross section has two factors, 
a function F(xt) {G{xt)) which 'scales', i.e. depends 
only on the ratio of momenta, and a dimensioned factor, 
l/Px )j where n equals 4 for QED, and for LO- 

QCD (Eq.0J, analogous to the l/g'* form of Rutherford 
Scattering. The structure and fragmentation functions 
are all in the F{xt) {G{xt)) term. The original high px 
measurements at CERN 28] and Fermilab [sj , showed 
beautiful xt scaling, but with a value of n = 8 instead 
of n = 4, for values of 3 < pt < 7 GeV/c. Later 
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measurements at larger px showed the correct scaling 
in agreement with pQCD and it was realized that the 
value n = 8 at lower values of pt and ^/s was produced 
by the fcrx smearing 0, 0|. More recently, the de- 
viation of tt'' and direct photon inclusive cross sections 
measurements from pQCD predictions has been used to 
derive the values of fcr required to bring the measured 
and smeared pQCD predictions into agreement. [26l| . 

A more direct method to determine fepy is to measure 
the acoplanarity of the dijet pair. Such measurements 
were originally performed at the CERN-ISR using two 
particle correlations Q, |M Hi 1^ ■ The same method 
will be used in the present work. 

Hard-scattering in p+p collisions at = 200 GeV is 
detected by triggering on a tt^ with transverse momen- 
tum pT^ > 3 GeV/c; and the properties of jets are mea- 
sured using the method of two particle correlations. The 
trigger tt*^ is a leading particle from a large transverse 
momentum jet while the associated particle comes from 
either the same jet or the away-side jet. We will ana- 
lyze an outgoing dijet pair, with trigger jet transverse 
momentum magnitude pTt which fragments to a trigger 
particle with transverse momentum pTt, and an away- 
side jet transverse momentum magnitude of pTa which 
fragments to a particle with transverse momentum ptb.- 
The average transverse momentum component of the 
away-side particle pVa perpendicular to trigger particle 
PTt in the azimuthal plane is labeled as Pout- If the mag- 
nitude of the jet transverse fragmentation momentum 
jr (Fig. ^a) is neglected, the magnitude of v^/cTy can 
be related to Pout- V^kxy = PoutpTa/PTa = Pout/za- 
Thus the measurement of Pout and the knowledge of the 
fragmentation variable (za) determines the magnitude 
of the parton's transverse momentum kx- 

The smearing of the steeply falling parton px spec- 
trum by the fcxx distribution tends to make the trigger 
jet transverse momentum pTt larger than the away jet 
transverse momentum pTa ■ The component of the net 
transverse momentum of the parton pair along the trig- 
ger direction is smeared by \/2fcTx such that: 

((?5Tt-PTaxf)=2(4x)=(4) • (5) 

For a flat px spectrum, the smearing would average to 
zero so that there would be no net shift in the transverse 
momentum spectrum: 

(pTt - Pt) = {pT - pTax) = . (6) 

However, due to the steeply falling px spectrum, the 
kxx smearing results in a net imbalance of the jet-pair 
towards the trigger direction. In the limit when kx is 
collinear with the trigger jet and with the requirement 
of the Lorentz invariance of s {Px=pTtpTa) it is easy to 
see that 

{pTt - Pt) = (^{Pt - pTa)) - ^(PTt - pTa) > . 
\PT / 2 

(7) 

We denote the imbalance of pxa and pxt by the quantity 

Xh = ipTa) / ipTt) ■ (8) 



Jet fragments have a momentum jx perpendicular to 
the partonic transverse momentum (Fig.^b). This vec- 
tor is again a two-dimensional vector with one compo- 
nent perpendicular to the jet transverse axis, px in the 
transverse plane and the other component perpendicu- 
lar to the jet transverse axis in the longitudinal plane 
(defined by the beam and jet axes). The component 
of jx projected onto the azimuthal plane is labeled as 
jxy The magnitude of (jxy), the mean value of jx 
projected into the plane perpendicular to the jet thrust 
(see App IA measured at lower energies has been 
found to be px independent and « 400 MeV/c, consis- 
tent with measurements in e~^e~ collisions [3j,|33j. 

This analysis uses two-particle azimuthal correlation 
functions to measure the average relative angles be- 
tween a trigger 7r° and an associated charged hadron. 
The angular width of the near- and away-side peak in 
the correlation function is used to extract the value of 
(jx) and x^^ (zt) a/ {k'^). An analysis of the associ- 
ated yields is used to extract the fragmentation func- 
tion which provides the (zt) and (za) values used for 
(k^) extraction. The details on the PHENIX experi- 
ment relevant to this analysis follow. 

III. EXPERIMENTAL DETAILS 

The PHENIX experiment consists of four spectrom- 
eter arms - two around mid-rapidity (the central arms) 
and two at forward rapidity (the muon arms) - along 
with a set of global detectors. The layout of the 
PHENIX experiment during RHIC Run-3 is shown in 
Fig.H 

Each central arm covers the pseudorapidity range 
I77I < 0.35 and 90 degrees in azimuthal angle (j). In each 
of the central arms, charged particles are tracked by a 
drift chamber (DC) positioned from 2.0 to 2.4m radially 
outward from the beam axis and 2 or 3 layers of pixel 
pad chambers (PCI, PC2, PC3 located at 2.4m, 4.2m, 
5m in the radial direction, respectively). Particle iden- 
tification is provided by ring imaging Cerenkov coun- 
ters (RICH), a time of flight scintillator wall (TOF), 
and two types of electromagnetic calorimeters (EM- 
Cal), lead scintillator (PbSc) and lead glass (PbGl). 
The magnetic field for the central arm spectrometers 
is axially symmetric around the beam axis. Its com- 
ponent parallel to the beam axis has an approximately 
Gaussian dependence on the radial distance from the 
beam axis, dropping from 0.48 T at the center to 0.096 
T (0.048 T) at the inner (outer) radius of the DC. A 
pair of Zero-Degree Calorimeters (ZDC) and a pair of 
Beam-Beam Counters (BBC) were used for global event 
characterization. Further details about the d esig n and 
performance of PHENIX can be found in |3 [s^l [s^ . 

A p + p data sample corresponding to an integrated 
luminosity 0.35 pb^^ at y/s — 200 GeV has been used 
for the present analysis. It contains a minimum bias 
(MB) sample of 121M events and a high-p^ triggered 
sample of 50M events. The MB trigger is obtained 
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FIG. 2: (color online) The PHENIX experimental layout for 
the Au+Au run in 2003. The top panel shows the PHENIX 
central arm spectrometers viewed along the beam axis. The 
bottom panel shows a side view of the PHENIX muon arm 
spectrometers and the position of the global detectors (BBC 
and ZDC). 



from the charge multiplicity in the two BBCs situ- 
ated at large pseudo-rapidity (ry « ±(3.0 — 3.9)). The 
BBCs were also used to determine the collision ver- 
tex, which is limited to a ±30cm range in this analysis. 
The high-pT trigger requests an additional discrimina- 
tion on sums of the analog signals from non-overlapping, 
2x2 groups of adjacent EMCal towers situated at mid- 
rapidity {\r]\ < 0.35) equivalent to an energy deposition 
of 750 MeV ^37] ■ The analysis has been performed sepa- 
rately on the two data sets and no trigger selection bias 
was found within the quoted errors. 

Neutral pions, which are used as trigger parti- 
cles, are detected by the reconstruction of their 77 
decay channel. Photons are detected in the EM- 
Cal, which has a timing resolution of « 100 ps 
(PbSc) and « 300 ps (PbGl) and energy reso- 
lution of aE/E=1.9%®8.2%/y/E{GeV) (PbSc) and 
aE/E^0.8%(B8A%/y/E{GeV) (PbGl). In order to im- 
prove the signal/background ratio we require the min- 
imum hit energy > 0.3 GeV, a shower profile cut as 
described in 38] , and no accompanying hit in the RICH 
detector, which serves as a veto for conversion electrons. 
A sample of the invariant mass distribution of photon 
pairs detected in the EMCal is shown in Fig. |21 
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FIG. 3; The measured 77 invariant mass distribution for 
pair pt in 4 < pr-y-y < 5 GeV/c. The peak is fitted with 
a Gaussian. The signal/background ratio within 2a of the 
mean ranges from ~ 6 at pt of 3 GeV/c up to « 15 at 
8 GeV/c. 



Charged particles are reconstructed in each PHENIX 
central arm using a drift chamber, followed by two lay- 
ers of multiwire proportional chambers with pad read- 
out Is^l • Particle momenta are measured with a resolu- 
tion Sp/p = 0.7% ® l.l%p (GeV/c). A confirmation hit 
is required in PC3. We also require that no signal in the 
RICH detector is associated with these tracks. These 
requirements eliminate charged particles which do not 
originate from the event vertex, such as beam albedo 
and weak decays, as well as conversion electrons. 

High momentum charged pions (above the RICH 
Cerenkov threshold) are identified using the RICH and 
EMCal detectors. Candidate tracks must be associated 
with a hit in the RICH [^^i which corresponds to a 
minimum momentum of 18 MeV/c for electrons, 3.5 
GcV/c for muons, and 4.9 GeV/c for charged pions. In 
a previous PHENIX publication we have shown 
that charged particles with reconstructed px above 4.9 
GeV/c, which have an associated hit in the RICH, 
are dominantly charged pions and background electrons 
from photon conversions albedo. The efficiency for de- 
tecting charged pions rises quickly past 4.9 GeV/c, 
reaching an efficiency of > 90% a.t px > Q GeV/c. 

To reject the electron background in the charged pion 
candidates, the shower information at the EMCal is 
used. Since most of the background electrons are gen- 
uine low pt particles that were mis-reconstructed as 
high Pt particles, simply requiring a large deposition of 
shower energy in the EMCal is effective in suppressing 
the electron background. In this analysis, a momentum- 
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FIG. 4: (color online) The probability distribution for 
charged pion candidates and electrons derived from the EM 
shower shape using identified electrons and pions. The in- 
tegrals have been normalized to one. 




15 

PT(GeV/c) 



FIG. 5: (color online) The raw charged pion transverse mo- 
mentum spectrum, with the final cuts applied. The level of 
the remaining background is estimated from an extrapola- 
tion from low-pT and is shown as a black line. 



dependent energy cut on the EMCal is applied 

£; > 0.3 + 0.15pT- (9) 

In addition to this energy cut, the shower shape in- 
formation js^ is used to further separate the broad 
hadronic showers from the narrow electromagnetic 
showers and hence reduce the conversion backgrounds. 
The difference of the EM shower and hadronic shower 
is typically characterized by a variable [s^ , 



E 



(10) 



where Ef^'^"-" is the energy measured at tower i and 
j^pred ^-j^g predicted energy for an electromagnetic 
particle of total energy E^'^"''' . 

In this analysis we use the probability calculated from 
this value for an EM shower, ranging from to 1 with 
a flat distribution expected for an EM shower, and a 
peak around for an hadronic shower. 

Figure 0] shows the probability distribution for pion 
and electron candidates, each normalized to one. The 
pion candidates were required to pass the energy cut 
and the electrons were selected using particle ID cuts 
similar to that used in The electron distribution 

is relatively flat, while the charged pion distribution 
peaks at 0. A cut of shower shape probability < 0.2 
selects pions above the energy cut with an efficiency of 
> 80%. Detailed knowledge of the pion efhciency is 
not necessary, since we present in this paper the per- 
trigger pion conditional-yield distributions, for which 
this efficiency cancels out. 

Since the energy and shower shape cuts are indepen- 
dent of each other, we can fix one cut and then vary the 
second to check the remaining background level from 
conversions. The energy cut in Eq. (^J is chosen such 
that the raw pion yield is found to be insensitive to 
the variation in the shower shape probability. Figure 



shows the raw pion spectra for EMCal-RICH triggered 
events as a function of pt, with the above cuts applied. 
The pion turn on from 4.9 — 7 GeV/c is clearly visible. 
Below Pt of 5 GeV/c, the remaining background comes 
mainly from the random association of charged parti- 
cles with hits in the RICH detector. The background 
level is less than 5% from 5—16 GeV/c, which is the pr 
range for the charged pion data presented in this paper. 



IV. CORRELATION FUNCTION 

The analysis uses two-particle azimuthal correla- 
tion functions between a neutral pion and an asso- 
ciated charged hadron to measure the distribution of 
the azimuthal angle difference Ac/) — <j)t — ipa (see 
Fig. O. Whenever a tt" was found in an event, the 
real, diVuncorr/rfA^, and mixed, dN^i^/ d/S^cj), distribu- 
tions for given pTt {t^^) and pTa. (charged hadron) were 
accumulated (left panel of Fig. EJ. Mixed events were 
obtained by randomly selecting each member of a par- 
ticle pair from different events having similar vertex po- 
sition. Then the mixed event distribution was used to 
correct the correlation function for effects of the limited 
PHENIX azimuthal acceptance and for the detection ef- 
ficiency, to the extent that it remains constant over the 
data sample. 

We fit the raw dA^uncorr/rfAc/) distribution with the 
product 



dNu 



dA(l) 
1 dN„ 



(11) 



TV dAcj) 



(Co + Ci ■ /ncar(A<^) + C2 • /away(A</')) 



where the mixed event distribution is normalized to 27r 
(TV = Y^dNl^^^JdAcf) see blue dashed line on the left 
panel of Fig. O, Co_2 are constant factors to be de- 
termined from the fit, /near(A(^) and /away(A(/)) are the 
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FIG. 6: (color online) An example of the correlation func- 
tions for 3< pTt <3.5 GeV/c and associated particles in 
1.4< pTa <5 GeV/c. (upper) Unnormalized pair-yield dis- 
tribution plotted with the fit function which is two Gaus- 
sians modulated by the pair detection efficiency derived from 
the mixed distribution (blue dashed line), (lower) Per vr'' 
trigger yield distribution corrected for the pair detection ef- 
ficiency. Dashed line represents the constant term in the 
fit. 



near- and away-side peak fit functions respectively. Tra- 
ditionally, the Gaussian functions, around = and 
around Ai^ = tt, are used for /near(A<?!)) and /away (A0). 
This leaves a total of five free parameters to be deter- 
mined - the areas and widths of the above two Gaus- 
sians: Yn; cn for the near-angle component and Yp, a a 
for the away-angle component and the constant term 
describing an uncorrelated distribution of particle pairs 
which are not associated with jets. However, the as- 
sumption of the Gaussian shape of the angular corre- 
lation induced by jet fragmentation is justified only in 
the high-pT region where the relative angles are small. 
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FIG. 7: (upper) Inclusive charged hadron efficiency correc- 
tion function, (lower) rj acceptance correction factor for loss 
of jet pairs outside the limited 7)-acceptance of the PHENIX 
experiment. 



In order to access also a lower px region we used an al- 
ternative parameterization of /ncar(A0) and /away (A0) 
which will be discussed later in the text. 

The normalized correlation function was constructed 
as a ratio of real and mixed distributions multiplied 
by 77-acceptance correction factor Ratj, divided by pr- 
dependent efficiency correction e{pT) (see left panel of 
Fig. 0) and divided by the number of tt" triggers. 

1 dN ^ RArj djVuncorr(A(/))/dA0 

NtriggdAcj) iVtrigge(pT) dNn,i^{A(f>)/dA(t) 

(12) 

The Rati correction factor which accounts for lim- 
ited 7] acceptance of the PHENIX experiment (see right 
panel of Fig. [TJ for the the near-side yield, with an as- 
sumption that the angular jet width is the same in A77 
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FIG. 8: (left) Measured yield of charged hadrons with away-side transverse momentum 1.4 < pTa < 5.0 GeV/c associated 
with a trigger tt" of transverse momenta given in in Table |I] (right) Measured yield of charged hadrons associated with a 
trigger tt" of fixed transverse momentum 3.0 < prt < 10.0 GeV/c and the away-side transverse momenta given in Table|I] 
The dashed lines corresponds to the fit of two Gaussian functions representing the trigger (t) jet and away-side (a) jet 
correlation. The x^(DOF) ctn and values extracted from these fits are tabulated in Table U 



TABLE I: The x^(DOF) ctn and values extracted for the correlation function shown in Fig.|Hl All units in rad and 

GeV/c. Only the statistical errors are shown. 





1.4 < pTa 


< 5.0 GeV/c 








3.0 < PTt < 10.0 GeV/c 




PTt GeV/c 


X=^(DOF=34) 




V (Pout) 


PTa 




x'(DOF=34) 


(Tn 


V(Pout) 


(a) 2.5 - 3.0 


69.4 


0.26 ± 4E-03 


1.17 ± 0.07 


(a) 1.0 - 


2.0 


188.4 


0.29 ± 4E-03 


0.87 ± 0.03 


(b) 3.0-3.5 


79.6 


0.24 ± 4E-03 


1.19 ± 0.05 


(b) 2.0- 


3.0 


63.2 


0.21 ± 3E-03 


1.16 ± 0.04 


(c) 3.5 - 4.5 


61.2 


0.22 ± 3E-03 


1.04 ± 0.04 


(c) 3.0 - 


4.0 


50.3 


0.16 ± 4E-03 


1.36 ± 0.06 


(d) 4.5 - 5.5 


52.7 


0.22 ± 6E-03 


1.08 ± 0.06 


(d) 4.0 - 


5.0 


63.2 


0.14 ± 4E-03 


1.69 ± 0.13 


(e) 5.5 - 6.5 


38.4 


0.20 ± 8E-03 


0.90 ± 0.06 












(f) 6.5 - 8.0 


31.6 


0.16 ± lE-02 


0.64 ± 0.06 













and in A0, can be written as 

1 



rO.7 



Arj 



acc{Ari)dr] 



(13) 



where acc{A'i]) represent the PHENIX pair acceptance 
function in IA77I. It can be obtained by convolving two 
flat distributions in | Aryj <0.35, so acc{AT]) has a simple 
triangular shape: acc{Ar]) — (0.7— |A?7|)/0.7. For the 
away-side yield the corresponding Rati is 

2(0.7) 



RAtj 



1-0.7 acc{Ari)d-q 



(14) 



^At; equals 2, because the pair efficiency has a trian- 
gular shape in \A'q\ <0.7, which results in 50% average 
efficiency when the real jet pair distribution is flat in 
IA77I <0.7. Normalized correlation functions for various 
PTt and PTa are shown in Fig. |S1 

For two particles with transverse momenta PTt^ PTa 
from the same jet, the width of the near-side correlation 
distribution can be related to the RMS value of the 
jx vector component, jTy, perpendicular to the parton 
momentum as 



4 = = 



yPTa.) \PTt, 



(15) 



10 



where we assume (j|.y)^p^^ and p^^ and thus the arc- 
sine function can be approximated by its argument and 
we can solve for ^ 



72. 



PTt PTa 



(16) 



In order to extract (|A:Ty|), or {k'^)^ we start with 
the relation between the magnitude of Pout,(see 

Fig.Hl) 



Pout = PTa sin A(/), 



(17) 



which is the transverse momentum component of the 
away-side particle pVa perpendicular to trigger particle 
PTt in the azimuthal plane (see Fig. and kTy'- 



where 



2(|fcTy|)V(bTy|)' +(bTy|>', (18) 



PTt ■ PTa. 
pit 



PTa COS ZaPTa 



PTt 



ZtPTt 



(19) 



represents the fragmentation variable of the away-side 
jet. 2, 3J We note however, that |3| explicitly neglected 
(zt) = (pTt/PTt) in the formula at ISR energies, where 
(zt) ~ 0.85, while it is not neghgible at -^5=200 GeV. 
Furthermore, as mentioned earlier, the average values 
of trigger and associated jet momenta are generally not 
the same. There is a systematic momentum imbalance 
due to /cT-smearing of the steeply falling parton momen- 
tum distribution. The event sample with a condition of 
PTt>PTa is dominated by configurations where the kr- 
vector is parallel to the trigger jet and antiparallel to 
the associated jet and (pTt — PTa) 7^ 0. Here we intro- 
duce the hadronic variable Xh in analogy to the partonic 
variable Xh 



PTa 
PTt 



Xh = Xh{{kl) ,Xh) 



{Pl 



'Tt/ 



(20) 



The detailed discussion on the magnitude of this im- 
balance is given later. In order to derive the relation be- 
tween the magnitude of Pout and kx let us first consider 
the simple case where we have neglected both trigger 
and associated (jx) (see panel (a) on Fig. In this 
case one can see that 



(IPoutI) |jrt=jra=0 = (boutl)oo — 

= V2{\kTy\) 



PTa, 



Ta) 



V2{\kTy\){zt)^ 



^ For relations between yJ{X'^) and (\Xy\), where X is any 2D 
quantity, see Add I A l1 



Rewriting the formula for pout in terms of RMS we get 

■ Xh 

Xh 



Ttj 



where we have taken (k1)~ 

However, the jet fragments are produced with finite 
jet transverse momentum jt- The situation when the 
trigger particle is produced with jTty> GeV/c and 
the associated particle with jTay^O GeV/c is shown in 
Fig-Ch- The Pout vector picks up an additional compo- 
nent 



(Pout) ljTt>0,iTa=0 



(Pout) 



00 



P\t 



(PTa ~ (Pout)oo) 



PTt 



(^Ity) 



P\t 



With an assumption of j^ty ^ PTt we found that 



(Pout) liTt>0, jra=0 — x\ 



1 



{ztY (4) + O'Tty) 



We include jTa in the same approximation, j^ay ^ PTa, 
i.e. collinearity of jTa and Pout with result 



(Pout) 



1 



(^ty(4)^ + (j|ty) 



+ OTay) (21) 



and we solve for (zt) \J (fc^) /^h 



Xh 



out) 



JTa 



^h WTtv 



If we assume no difference between j^t and j'^a then we 
have 



(zt(fcT,Xh)) V(fc 

Xh(A;T,Xh) 



1 



Xh 



y(pLt)-(j|y)(l+:^g) 



(22) 

All quantities on the right-hand side of Eq. (|22|l can be 
directly extracted from the correlation function. The 
correlation functions are measured in the variable Ac/i 
in bins of pTt and pTa (e.5. see Fig. ISjl, and the rms of 
the near and away peaks ctn and cta are extracted. We 
tabulated ctn and cta for many combinations of p^t and 
PTa (see Fig. and Fig. UHl . 

Initially, we used the approximation VvPout) ^ 
PTa, sin a A in Eg. 1221 However, we have noticed that this 
approximation and other approximations for \/ (Pout) 
proposed e.g. in reference |43| (see appendix lA 2|l are 
inadequate for cta > 0.4 radians, so we don't use cta to 
calculate kx- 

We extract y/ (pl^t) directly for all values of pTt PTa 
(even for wide bins in pTa using the (pTa) of the bin) by 
fitting the correlation function in the tt/2 < A(f> < 3tt/2 
region by 



dAcj) 



3-!t/2 
Tr/2 



_ dN dpout 
dpout dAcf) 
PTa cos A0 



(23) 



VMSjErf(-^ 



■ exp 



Ptjj sin^ A(f)^ 

2{pLt) , 
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TABLE II: Measured widths of the near- and away-angie tt" — correlation peaks for various trigger particle momenta 
Only the statistical errors are shown. 



PTt=3.39 


GeV/c 


PTt=4.40 


GeV/c 


PTt=5.41 


GeV/c 


PTt=6.40 


GeV/c 


PTa CTN rad 


(JA rad 


PTa o-N rad 


a A rad 


PTa o-N rad 


(TA rad 


PTa O-N rad 


OA rad 



1.59 0.27 ± 0.01 0.58 ± 0.05 1.72 0.28 ± 0.02 0.50 ± 0.03 1.51 0.26 ± 0.01 0.49 ± 0.03 1.34 0.40 ± 0.03 0.68 ± 0.05 

1.84 0.24 ± 0.01 0.52 ± 0.03 2.14 0.18 ± 0.01 0.47 ± 0.06 2.22 0.21 ± 0.02 0.39 ± 0.05 1.64 0.30 ± 0.02 0.58 ± 0.05 

2.22 0.23 ± 0.01 0.52 ± 0.03 2.53 0.20 ± 0.01 0.47 ± 0.04 2.88 0.17 ± 0.01 0.37 ± 0.05 1.94 0.23 ± 0.02 0.52 ± 0.06 

2.73 0.19 ± 0.01 0.46 ± 0.04 3.17 0.16 ± 0.01 0.38 ± 0.04 4.01 0.14 ± 0.02 0.34 ± 0.07 2.29 0.23 ± 0.02 0.40 ± 0.03 

3.24 0.19 ± 0.01 0.47 ± 0.04 4.36 0.14 ± 0.01 0.39 ± 0.07 2.74 0.17 ± 0.01 0.41 ± 0.05 

3.93 0.17 ± 0.01 0.41 ± 0.03 3.36 0.17 ± 0.02 0.36 ± 0.04 

5.04 0.12 ± 0.01 0.38 ± 0.05 



where we assumed a Gaussian distribution in Pout- 
We still use a Gaussian function in Acj) in the near angle 
peak to extract ^ {Jt)- The a/ {pint) values extracted 
from the fit of the functional form (|23|l are shown in 
Fig. El and Fig.ini 
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FIG. 9: (color online) (top) The width of the near-side peak 
(7N with PTa for various values of pTt as indicated in legend, 
(bottom) The width of the far-side peak ax with pTa for the 
same prt selection. 

The per-trigger yields as a function of prt for fixed 
associated 1.4 < pTa, < 5.0 GeV/c bin are shown in 
Fig. ^1 There is a distinct behavior of the near-side 
yield which varies with trigger p^t much less than the 
away-side yield. For the away-side, this reflects the fact 
that the particle detected in the fixed associated bin 
are produced from the lower z region of the fragmenta- 
tion function for events with higher trigger pxt ■ For the 
near-side jet, this multiplicity increase is reduced due 
to the fact that with increasing pxt the near-side jet 
energy increases; however, at the same time the larger 
fraction of this energy is taken away by the more ener- 
getic trigger particle. Thus the relative change in z is 
smaller on the near-side. 
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FIG. 10: (color online) The near-side (squares) and away- 
side (circles) width as a function or trigger- tt" prt- The asso- 
ciated charged particle momenta are in the 1.4 < pTa < 5.0 
GeV/c region. The curves are from a PYTHIA calculation 
with the values of kx indicated. The data values are given 
in Table ITTTl 



TABLE III: The ctn and a a values shown in Fig. [TOl All 
units in rad and GeV/c. Only the statistical errors are 
shown. 



PTt 



(TN 



UA 



2.23 
2.72 
3.22 
3.89 
4.90 
5.91 
7.23 



0.247 
0.227 
0.235 
0.215 
0.210 
0.197 
0.185 



0.002 
0.003 
0.004 
0.004 
0.006 
0.009 
0.012 



0.565 
0.548 
0.521 
0.464 
0.431 
0.396 
0.350 



0.013 
0.014 
0.016 
0.014 
0.020 
0.025 
0.028 



In order to extract (zt) and Xh knowledge of the frag- 
mentation function is needed; a detailed discussion is 
given in following sections. 
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FIG. 11: Extracted values of ^/(J^ for 3.0 < prt < 4.0 
and 5.0 < prt < 10.0 GeV/c for various values of pTa using 
the direct pout extraction method based on fitting the away- 
side peak by Eq. 12311 . 
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FIG. 12: (solid circles) Extracted values of (Paut) fo'' 
1.4 < pTa < 5.0 GeV/c for various values of prt from 
Eq. 1231 . (open diamonds) Indirect pout=PTaSin((7A) values. 



(zt) ^ (fcy) results 

The measurement is performed in two different kine- 
matical regimes; first the transverse momentum of the 
trigger particle, pTt, is fixed and the peak width is mea- 
sured for different values of associated particle trans- 
verse momenta pTa (Fig- 19 • (Note that in the region 
of overlap, the data are in excellent agreement with a 



TABLE IV: 
Fig. m All 
are shown. 



The ^ (Pout) values shown in Fig. and 
units in GeV/c. Only the statistical errors 



1.4 

PTt 



< PT, 



< 5.0 



P out 



3 < PTt < 4 



(Pout 



5 < PTt < 10 

PTa y( 



;Pout) 



2.23 
2.72 
3.22 
3.89 
4.90 
5.91 
7.24 



1.315 
1.250 
1.182 
1.011 
0.953 
0.868 
0.798 



0.043 
0.046 
0.049 
0.038 
0.052 
0.064 
0.068 



1.72 0.996 ± 0.056 1.85 0.960 ± 0.102 



2.22 1.244 ± 
2.73 1.222 ± 

3.23 1.496 ± 
3.93 1.793 ± 
5.04 1.675 ± 



0.079 
0.095 
0.105 
0.115 
0.141 



2.24 
2.73 
3.44 
4.65 



1.100 ± 0.103 
1.088 ± 0.110 
1.285 ± 0.136 
1.268 ± 0.210 



TABLE V: The near and away side conditional yield per 
number of triggers for 1.4 < pTa < 5.0 GeV/c shown in 
Fig. 1131 All units in rad and GeV/c. Only the statistical 
errors are shown. 



PTt 



Ya 



2.23 
2.72 
3.22 
3.89 
4.90 
5.91 
7.24 



1.911 
1.863 
2.032 
1.966 
2.120 
2.153 
2.174 



0.018 
0.022 
0.032 
0.033 
0.061 
0.098 
0.125 



1.717 
1.908 
2.130 
2.360 
2.611 
2.992 
3.690 



0.044 
0.055 
0.071 
0.074 
0.123 
0.196 
0.242 
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FIG. 13: Measured yield of charged hadrons associated with 
one trigger tt" with transverse momenta indicated in TableQ] 
and associated charged hadron with 1.4 < pt<5.0 GeV/c. 
Dashed lines represent the linear fit. 
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TABLE VI: The ^/{jf) values shown in Fig.[T31and Fig.lTKl 

All units in rad and GeV/c. Only the statistical errors are 
shown. 





1.4 < PTa < 5.0 




3 < PTt < 4 


PTt 




PTa 




3.22 


587 ± 009 


1.72 


0.562 ± 0.011 


3.89 


0.577 ± 0.009 


2.22 


0.597 ± 0.014 


4.90 


600 ± 017 


2.73 


572 ± 017 


5.91 


0.596 ± 0.026 


3.23 


0.590 ± 0.020 


7.24 


0.597 ± 0.038 


3.93 


0.603 ± 0.017 


8.34 


0.632 ± 0.085 


5.04 


0.506 ± 0.029 




4 < pTt < 5 




5 < PTt < 6 


PTa 




PTa 




1.72 


0.643 ± 0.036 


1.52 


0.529 ± 0.030 


2.14 


0.492 ± 0.032 


2.22 


0.581 ± 0.049 


2.53 


0.608 ± 0.035 


2.88 


0.590 ± 0.047 


3.17 


0.590 ± 0.032 


4.01 


0.603 ± 0.063 


4.36 


0.631 ± 0.052 
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FIG. 14: (color online) (jt) values calculated according 



Eq. I|16|l . The dashed line represents the 0*' 
mial fit in the 1.5 < pTa < 5 GeV/c region. 



-order polyno- 



previous measurement ji^.) In the second case, parti- 
cle pairs with a fixed associated bin 1.4 < pTn < 5.0 
GeV/c and various prt are selected fFig. llO|l . It is evi- 
dent that both near and away-side correlation peaks in 
all cases reveal a decreasing trend with pTa and prt ■ 

However, the asymptotic behavior of ctn and (Ta is 
different. Whereas the magnitude of ctn, according to 
Eq. Hlt)|l . should vanish for large values of pTt and PTa, 
the (7A according to Eq. should be approximately 
constant around x^^ (zt) \/ (fcf.) /prt for large values of 
PTa- The (zt) and ih quantities are implicitly pxa de- 
pendent, however, their ratio is roughly ~ 1 so that the 
asymptotic value of ""Ajp^ ^y/Wr) /PTt- 

Extracted values of function of pTa accord- 

ing to Eq. (|16|l are shown in Fig. All values 
are constant in the explored region (pTa>l-5 GeV/c). 



FIG. 15: Averaged values of (jy) in (1.5< pTa <5 GeV/c) 
as a function of the trigger transverse momentum pTt (solid 
triangles). The CCOR values measured at ^/s=62.4 GeV 
shown by open triangles. 

It is expected that can not remain constant for ar- 
bitrarily small values of pxa because of the phase space 
limitation. In the region where PTa< VUr): the magni- 
tude of the jx-vector is truncated, similar to the "Seag- 
ull effect" j4^J . Since the \/Jj^ values are on the order 
of 600 MeV/c, we assume that the phase space limi- 
tation can be safely neglected for pTa>l-5 GeV/c and 
we extract the values of -y/ (j^) averaged over pTa and 
PTt (see Fig. HH) 



^JiJl) = 585 ± 6(stat) ± 15(sys) MeV/c (24) 

The systematic error originates from the finite mo- 
mentum resolution and Eq. I|16() where we assume that 
the arc-sine function can be approximated by its argu- 
ment. For the angular width of the near angle peak (see 
Fig. Inland Fig. I1U|) it corresponds to an uncertainty of 
order of 3%. 

The independence of {Jt} on either pTt or ^/s has 
been observed by the CCOR experiment in the range 
Vs=31-62.4 GeV P. The values at Vs=62.4 
GeV (open triangles on Fig. I15|) are systematically 
larger then values found in this analysis. The discrep- 
ancy should not be taken as significant, as CCOR used 
a slightly different technique than in this paper. CCOR 
extracted the {Jt} values from measurements of (|pout|) 
for different values of the x e variable Eq. H19() . Accord- 
ing to Eq. the (|pout|)^ magnitude should depend 
linearly on x'^; and the (|jTy|) value was extracted from 
the intercept of the (pout)^(2^£') fit at xe=0, rather than 
from a measurement of ctn. 

Knowing the andpout values, we used Eq. to 
determine ijj^^ (zt) \/{^^ (see Fig.lTCland Fig.[T7|). The 
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FIG. 16: (color online) iih ^{zt) ^ values calculated 
according to Eq. 1221 for trigger tt" in 3 < pTt < 4 GeV/c 
and 5 < prt < 6 GeV/c as a function of pTa. The system- 
atic uncertainties are indicated by the shaded boxes. 



systematic error was estimated with Monte Carlo simu- 
lations to be on the order of 5%. The main source of sys- 
tematic error originates from the assumption (Eqs. (|15|l 
and (|21|l 'l of the relative smallness of (ij-), coUinearity 
between pout and jTay and from the limited momentum 
resolution discussed in section lllll 



The PTa dependence of the extracted (zt) VV^r) 
(Fig. I16|) reveals a strikingly decreasing trend. It was 
originally expected that by fixing the value of pTt , the 
kinematics of the hard scattering (i.e. pTt—PTa) would 
be fixed, independently of the value of PTa.- Various 
values of pTa. would then sample the pTa. fragmentation 
function, and the value of x^^ (zt) \/(fcfT '^as expected 
to be constant. It is evident that this assumption is not 
quite correct. 

A similar line of argument applies also for the ris- 
ing trend when pTa is fixed and prt varies (Fig. [T7I) . 
It is interesting to note that the CCOR Vv^r) values 
measured at ^/s=62.4: GeV (open triangles on Fig. ll7|l 
reveal a similar rising trend. However, the rising trend 
of x^^ (zt) a/ {k'^) withpTt and falling withp^a suggests 
that the variation of V\MT with pTt seen by the CCOR 
collaboration 0] may be indicative of the {ztixj^^ varia- 
tion which was there neglected.^ In order to understand 
variation of (zt) and ih we have to explore the process 



^ Note, however that the method was different. CCOR deter- 
mined Jt and from the slope and intercept of Eg. I21l with 
respect to pxsu S't each value of , with the implicit assump- 
tion that {zt)xj^^ = 1. 



FIG. 17: (color online) The same calculation according 
Eq. ll^ for fixed associated bin 1.4 < pTa < 5.0 GeV/c 
as a function of prt- The systematic errors indicated by 
colored rectangles. The {k^) results obtained by CCOR 
collaboration at y^=62.4 GeV |l| displayed by solid trian- 
gles. 



of dijet fragmentation. 



V. FRAGMENTATION FUNCTIONS 

We have shown in Eq. (|22|l that the width of the away 
side correlation peak 

is related to the product of x^^ {zt) \/ (fcy)- In order 
to evaluate (zt), knowledge of the scattered parton p^ 
spectrum and fragmentation function is required. 

Fragmentation functions from e+e~ collisions, 
weighted by the appropriate hard-scattering constituent 
cross-sections and evolution could in principle be 
used. However, it was originally thought that the shape 
of the fragmentation function could be deduced from 
present measurements using the combined analysis of 
the inclusive trigger pxt and associated particle pxa. dis- 
tributions. Although this idea turned out to be incor- 
rect, we will follow this line of reasoning for a while as 
it is instructive. 

Generally, the invariant cross section for inclusive 
hadron production from jets can be parametrized in 
the following way. First, we assume that the number 
of parton fragments (consider only pions for simplicity) 
at a given px corresponds to the sum over all contribu- 
tions from parton momenta, px from pT < pr < \/s/2. 
The joint probability of detecting a pion with pt= zpx 
originating from a parton with px can be written as 



prdprdz prdpr 



X DUz) 



15 



(25) 



Here we use fqipr) to represent the final state scattered- 
parton invariant spectrum daq/pTdpT and D1{z) to 
represent the fragmentation function. The first term 
in Eq. (|25|l can be viewed as a probability of finding 
a parton with transverse momentum pT and the sec- 
ond term corresponds to the probability that the par- 
ton fragments into a particle of momentum pt— zpx- 
With a simple change of variables from px to pt= zpx, 
we obtain the joint probability of a pion with px which 
is a fragment with momentum fraction z from a parton 
with pt~Pt/z: 



PTdpTfiz 



(26) 



The Pt and z dependences do not factorize. However, 
the Pt spectrum may be found by integrating over all 
values of pt > Pt to pTmax = ^1^^ which corresponds 
to values of z from xt = '^■Pt / Vs to 1 . 



1 da,, 
Pt dpT 



(27) 



Alternatively, for any fixed value of pT one can evaluate 
the {z{pt)), integrated over the parton spectrum: 



{z{vt)) = 



jl^zDljz) f,{pT) {pt/z) 
jLDUz)fq{pT/z)^ 



dz 



(28) 



From the scaling properties of QCD and from the 
shape of the 7r° invariant cross section itself, which is a 
pure power law for pr > 3 GeV/c [s^l, one can deduce 
that fqipT) should have a power law shape, fq{pT) = 
Ap^". In this case the hadron spectrum also has a 
power law shape because 



1 dUTr 

Pt dpT 



ADl{z)-{ 



PT.-ndz 



Pt 



Dl{z)-z''-'dz 



(29) 



and the last integral depends only weakly on pT due 
to the small value of xt- For small parton pT (below 
3-4 GeV/c) the power law shape is no longer valid, 
but the region pT < 3 GeV/c is outside the scope of 
this paper. The fq{pT) should also diminish for very 
high pt^ where the phase space available for 

hard parton production diminishes, again not relevant 
for the present purposes. 

We used the power law parameterization for the final 
state scattered-parton invariant spectrum fq{pT) oc p^" 
where n is a free parameter which can be determined 
from the fit of Eq. (|27|l to the measured 7r° cross section. 
There is, however, one more missing piece of informa- 
tion - the shape of the fragmentation function D^. In 
an attempt to extract this information from the data, 
we have analyzed associated x^j-distributions, as shown 
in Table EUl 



A. 'Scaling' variable xe 

It was expected !^ that the xe variable, defined by 
Eq. I|19|) . to first order, approximates the fragmentation 
function in the limit of high values of pTt , where there 
is sufficient collincarity between the trigger particle and 
the fragmenting parton. In this case where jr^Prt and 
fcr^PTt one can assume that PTt=PTt / zt and x-e zt = 
ih PTa COS A(^/pTa — -^aj and thus the slopes oiD{za) 
and xe are related as 



{xe) {Zt) x^ - 



(30) 



The xe distributions of particles associated with trig- 
ger particles in the 3-8 GeV/c range of transverse mo- 
mentum are plotted in Fig. 1181 The dashed lines repre- 
sent exponential fits. The slopes of these exponentials 
range from —5.8 (3 < pTt < 4 GeV/c) to —7.8 (open 
symbols on Fig. lT^ . This is qualitatively and quantita- 
tively different from the similar measurement done by 
CCOR collaboration at -^5=62. 4 GeV where the slopes 
of exponential fits to the xe distributions were found 
to be w —5.3 and independent of the trigger transverse 
momenta. That observation also supported the hypoth- 
esis of the xe distribution being a good approximation 
of the fragmentation function. We also note that the 
xe distributions are not quite exponential and at large 
values of xe there is a tail similar to the power law tail 
of the single inclusive pt distribution. 
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FIG. 18: Tlie distribution of associated particles with xe 
variable for various trigger particle prt indicated in the leg- 
end. Exponential fits indicated by dashed lines. 

The reason why the xe distributions do not have the 
same slope for different pTt and why there is a "power 
law" tail at large xe is the same as that which causes 

^ (-^t) \/ {k"^) to decrease with the associated particle 
transverse momentum. It turns out that by sampling 
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TABLE VII: Measured xe distributions associated with various transverse momenta of the trigger tt . Only the statistical 
errors are shown. 

PTt=3.39 GeV/c ^^^=4.40 GeV/c prt^dAl GeV/c pTt=6.40 GeV/c pTt=7.39 GeV/c 
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different regions of pTa for fixed prt, the average mo- 
mentum of the parton fragmenting into a trigger parti- 
cle, (zt), also changes. This kind of trigger bias causes 
the hard scattering kinematics, the value of py, to not 
be fixed for the case where prt is fixed but pxa, varies. 
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FIG. 20: The distribution of associated particles with 
PTa/pTt variable for various trigger particle prt indicated in 
the legend. The distribution were fitted in the limited range 
0.2 < pTm/pTt < 0.4 by an exponential function (dashed 
lines) . 



FIG. 19: The negative slope parameters extracted from the 
fit of a plain exponential function into a xe (see Fig. 1181 and 
PTm/pTt (see Fig. I20II distributions. 

Taking this into account, one can not treat the asso- 
ciated Xe distribution as a rescaled fragmentation func- 
tion, but rather as a folding of the two fragmentation 
processes of trigger and associated jets. The same line of 
arguments applies also for other two-particle variables, 
e.g.pTa,/PTt, 45] used for an approximation of the frag- 
mentation variable z (see Fig. l20|) . The negative slopes 
of an exponential fit in the 0.2 < pTa/prt < 0.4 range 
(solid symbols on Fig. I19|l are, within the error bars. 



the same as for x^- 

In conclusion: the slope parameters extracted from 
associated x-^ distributions reveal the rising trend with 
PTt which reflects the fact, that the different pxa, sam- 
ples not only different but also different z^. 

The description of an associated distribution detected 
under the condition of the existence of a trigger particle 
requires an extension of the formulae discussed inlvland 
is a subject of the next section. 
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VI. DUET FRAGMENTATION 



For the description of the detection of a single particle 
which is the result of jet fragmentation, recall Eq. (|25|) 



dpxdzt 



(31) 



where we have now explicitly labeled the z of the trigger 
particle as zt, and defined 



^qiPr) = PT fqiPr) = ^ 
dpT 



(32) 



When kx smearing is introduced, configurations for 
which the high pT parton pair is on the average mov- 
ing towards the trigger particle are favored due to the 
steeply falling px spectrum, such that: 

{pTt -j5t) ^ ^ {pTt -PTa) = s(^t) 

with small variance a^, and we explicitly introduced 
PTt and pth to represent the transverse momenta of 
the trigger and away partons. The single inclusive pxt 
spectrum is now given by 



dpTtdzt 



I];(pTt) X Dlizt) 



(33) 



where the trigger parton pxt spectrum after fcy smear- 
ing is 



KiPTt) =PTt f'nipTt) 



dog 
dp 



Tt 



(34) 



Then, the conditional probability for finding the away 
side parton with pTa and z^, given prt (and zt), is: 



dP{pT^, Za 



dpTs.dz^ 



C{pTa.jPTt, kT)Dl{z^ 



PTt 



where C(pTs.,PTt, kx) represents the distribution of the 
transverse momentum of the away parton pTa, given 
pTt and kr, which can be written as: 



CipTa,PTt, fcr) 
1 



(35) 



exp 



[PTa 



Then 



2a2 



dP{pTa, Za) 



dpTtdztdpTadzi, dpTtdzt 



dpTadz^ 



In general, as/s{kT) is small (see section lVI Bp so that 
C{pTa.iPTt: kr) is well approximated by a (5 function and 
we may take 



PTa = PTt - 2s{kT) = XhPTt 



so that 



where 



dpTtdztdza, 



Y.'g{pTt)DUzt)Dl{z,) 



PTa PTa 



ZtPTa 



PTa XY,PTt XhPTt 

Changing variables from prt, zt to pTt, zt as above, and 
similarly from Za to pTa, we obtain 



1 ^, rPTt 



dpTtdztdpTa XhPTt ' Zt 



Y.'g{^-±l)Dl{zt)Dl{^) (36) 



ZtPTa 



where for integrating over zt or finding (zt) for fixed pTt, 
PTa, the minimum value of zt is z""" = IpTtj \fs = a^Tt 
and the maximum value is: 



max = j-i^PZl 



PTa 



2^ 



where a:h(pTt,PTa) is also a function of kx (Eq. 1(201) )• 

Thus, in order to evaluate a;h(pTt,PTa) for use in 
Eq. H3t)|l . fcy must be known. We attack this prob- 
lem by successive approximations. First we solve for 
kT and D%(t.) assuming = 1 as done at the ISR 
where the smearing correction was small. Then we 
solve for x\Sj>T\.-,PTa) with this value of fc^ and iterate. 
On the first solution we solve only for Y!g{pTt) while 
on the iteration we include the fcr smearing to solve 
for the unsmeared parton spectrum Eq(pT)= PT fqipT) 
(Eq. (|22J). 



A. Sensitivity of the associated spectra to the 
fragmentation function 

As discussed in section IVAI the associated xe distri- 
bution was thought to approximate the fragmentation 
function of the away jet. Equation 1361 can be trans- 
formed to the X E distribution at fixed pxt with a change 
of variables from pTa to Xe followed by integration over 

Zt- 



d^a 



dpTa 



d^a 



dpTtdxE dxE dpTtdpTa 

PTt 

— I Dl{zt)Dl{- )I] ( )dzt. 

Xh JxTt XY,PTt Zt 



(37) 



We at first attempted to solve for the fragmentation 
function by simultaneous fits of the measured x e distri- 
butions to Eq. H37|l constrained by a fit of the inclusive 
invariant tt" cross section to Eq. (|27|) . There were diffi- 
culties with convergence. 

The reason for the lack of convergence became appar- 
ent when we calculated xe distributions according to 
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FIG. 21: (color online) The same xe distributions as on 
Fig. 1181 shown with calculations according to Eq. I|37|l for 
quark (solid lines) and for gluon (dashed lines) D{z). An ex- 
ponential approximation was used and the slopes for quark 
and gluon fr agm entation function were obtained by fitting 
to LEP data and ^ (see Fig.l^. 



Eq. 1)3 7|l (Fig. I21|l for two different fragmentation func- 
tions corresponding to quark and gluon jet fragmenta- 
tion. A simple exponential parameterization was used 
and the slopes were obtained from the fit to the LEP 
data "i^, 'It'I (Fig. For quark and gluon jets, we 
found Dq{z) K, exp(— 8.2 • z) and Dg{z) k, exp(— 11.4 • z) 
respectively. For the parton final state spectrum, we 
used c>c p^^. It is evident that the x-e, distribu- 
tions calculated for the quite different quark and gluon 
fragmentation functions do not differ significantly (the 
difference between solid and dashed lines on Fig. I21() . 
Clearly, the xe distributions are rather insensitive to 
the fragmentation functions of the away jet in contradic- 
tion to the previous conventional wisdom. The evidence 
of this explicit counter example led to attempts to per- 
form the integrals of Eq. H3()|l and Eq. (|33|l analytically 
which straightforwardly confirmed the observation that 
the Xe distribution is not sensitive to the fragmentation 
function. 

If the smeared trigger parton spectrum is taken as a 
power law, 



S'(^)^A(^)-("-i) 



and the fragmentation function as an exponential, 
D{z) = B exp(— 6z), then the integral of Eq. (|36|l over 
Zt becomes: 



FIG. 22: Fragmentation function measured in e^e~ col- 
lisions at y^=180 GeV for gluon and quark jets. The 
solid and dashed lines represent the exponential fit in the 
0.2 < z < 1 region. 



^4rr^^ dztzr^eM-bzdi-^^)] 

which is an incomplete gamma function. Since ~ 1, 
we make the assumption that it is constant. Similarly, 
the integrals of Eqs.[2nilSSlare also incomplete gamma 
functions: 



dpTt 



I dztz" ^exp[-6zt] 



(39) 



A reasonable approximation for the inclusive single, 
and two particle cross sections is obtained by taking 
the lower limit to zero and the upper limit to infinity, 
leading to the replacement of the incomplete gamma 
functions by gamma functions, with the result that: 



r(n)S2 A 



1 



dpTtdpTa. 6" ih PTt (1 + if^ 
rfcr^ r{n - 1) AB 



PTt ' 



dp' 



Tt 



PTt 



(40) 



(41) 



where T{n) = {n - l)r{n - 1). 

The conditional probability is just the ratio of the 
joint probability Eq. (|in)l to the inclusive probability 
Eq. (||l}, or 



dP^ 



dpTa 



B{n-1) 1 



bpTt Xh (1 



In the coUinear limit, where pTa 



XhPTt ' 

^X^PTt 



(42) 
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The only dependence on the fragmentation function, 
in this approximation, is in the normahzation constant 
B/h which equals (m), the multiplicity in the away-jet 
from the integral of the fragmentation function. The 
dominant term in Eq. (|43() is the Hagedorn function 
1/(1 -l-XE/ih)", so that at fixed pxt the distribution 
is predominantly a function only of xe and thus does ex- 
hibit 'xe' scaling. Also, the Hagedorn function explains 
the "power law" tail at large xe noted in section IV^ 
The reason that the xe distribution is not very sensi- 
tive to the fragmentation function is that the integral 
over zt for fixed pTt and pTa, (Eq. (|2Hl) is actually an 
integral over the jet transverse momentum . How- 
ever since both the trigger and away jets are always 
roughly equal and opposite in transverse momentum, 
integrating over simultaneously integrates over pT^ , 
and thus also integrates over the away jet fragmentation 
function. This can be seen directly by the presence of 
Zt in both the same and away fragmentation functions 
in Eqs.|211and|23 so that the integral over zt integrates 
over both fragmentation functions simultaneously. 



falling parton distribution function and the fragmen- 
tation function is peaked at z « 1 preferring "small" 
parton momenta. We have developed more accurate de- 
scription of the conditional yields taking into account 
the kx smearing. 




FIG. 23: (color online) Back-to-back partons in hard scat- 
tering rest frame (blue dashed arrows) with four-momenta 
(Pt,0,0,pt) and (-pt,0,0,pt) in (-,-,-,+) metrics moving 
along pn (j3n = PTpair) for au event where detection of prt 
and PTa is required (the jr contribution is neglected). The 
Pt%> Pth condition implies that the events with pn pointing 
more in the direction of pTt are selected. 



B. 



feT smearing 



In order to evaluate XhipTt,PTa,) and fcy must be 
known. We attack this problem by successive approxi- 
mations: first we solve for kx assuming a;^ = 1 as done 
at the ISR, where the smearing correction was small. 
Then we iterate for finite kx- The Gaussian approxi- 
mation for the smearing function Eq. (|35|l does not work 
so well in the low px region. The product of the steeply 



Let us consider the configuration depicted on Fig. 1231 
The two back-to-back partons in s frame undergo the 
Lorentz boost determined by net pair momentum 

Pn = PTpair = Pxt + PTh = ^Tt + kxs^ (44) 

If we denote an angle between the unsmeared parton 
momentum and /cr-vector (or p„) as a (see Fig. ESI) then 
we can write the conditional probability distribution of 
trigger parton momenta, pxt, as 



dpTtdctdpT 



PTt ,PTa. 



PTt Pxt PT^{n) Jrt) 



r 



(45) 



where G[pn) — exp(— p^y2 (fc|,)) describes the Gaus- 
sian probability distribution of the net pair momentum 
magnitude distribution, T,q{px) is the unsmeared par- 
ton momentum distribution, D'^ is the fragmentation 
function and rt = ipTt, iPtPTi kx) is the phase space 
vector. The pxt is chosen to be an integration variable 
and pxa is fully determined by given values of pxt , Pt , 
angle </) and by the requirement of Lorentz invariance. 

In order to evaluate (zt(kx)) and 
Xh(kx)\ we have to evaluate first the par- 

ton distribution for events where given pxt and pxa. 
are detected. This conditional cross section can be 
expressed as a definite integral over the unobserved 
variables (j) and px (see Fig. ESI 



da 



dp 



'Tt 



= 2 



d^a 



dpxtdpxd<j) 



dcjydpj 



PTt PTt Jo 



v/i/2 



Sg(pT)x 



X / Urt)Gipr.in))-DU^ 
Jo 



PTa 



PTairt) PxAn) 



(46) 

The da/dpxal distribution can be derived from 



Eq. igni just by rotation pxt PTa and pxa. PTt- 
The (zt(kx)) and Xh(kx)\ quantities can 

then be evaluated as 

(-*(^-)) lpTt,PTa - m (47) 

where 

X /;pnG(p„(r.O) • d4> dpx dz. 
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and rzt = (prt/ zt, (I),Pt, kx). The ^h(fcT) 1^^^ ^^^ is 
evaluated as 



XhikT)\ 



t,PTa (py^) Xa{Q)Xt{l) 



(48) 



where 



X /;p,(rl) G(p,(rl)) . DUj^)^) d<t> dprPrt 

Xa{n) = ^;^frT:'Dl{f^Jjf'^^,{pT)y. 

X /;pn(ra) G(p„(r,)) • Dlij^)^;^^ dc^ dprPT. 
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FIG. 24: PYTHIA simulated average momentum unbalance 
for the associated particles in 3.0 < pTa < 4.0 GeV/c bin 
and calculated according Eq. 1481 . The two vertical dashed 
hne indicates the range where prt and pTa bins are equal and 
the parton momenta unbalance vanishes (fixed correlations). 

We have tested the above formulae on PYTHIA sim- 
ulation. We have generated events with ^/{^^=3 
GeV/c and evaluated the partons' momenta unbal- 
ance variation with pTt for fixed 3 < pxa < 4 GeV/c 
bin. The results from the PYTHIA simulation (solid 
point on Fig. I24|l are compared to calculation based on 
Eq. ll^ (solid line on Fig. |^ . The magnitude of mo- 
mentum unbalance saturates at pTt~10 GeV/c around 
V i^T^) ^^'^ then starts to decrease. The maximum 
value depends on the the kx magnitude and on the 
asymmetry between pTt and pTa- Eventually, the un- 
balance should vanish at high pxt as a consequence of 
T,q{pT) flattening. 

The comparison of (zt) and (za) found in PYTHIA 
and derived according to Eq. l|17|l is shown in Fig. [23 
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FIG. 25: Average z of a trigger and associated particle as 
a function of pTa from PYTHIA and according Eq. 1471 . 



The overall agreement between the PYTHIA simula- 
tions and the calculation is excellent. The small devia- 
tions may be attributed to the fact that in the PYTHIA 
simulation, 1 GeV/c-wide bins were used for trigger and 
associated particle identification, whereas the calcula- 
tion was performed for fixed values of pTt and pTs^■ 

The last missing piece of information needed before 
solving Eq. H22|l is the fragmentation function D1 and 
unsmeared T,q(pT). The description of how this knowl- 
edge was extracted from the data is a subject of next 
section. 



VII. CORRECTED (kr) RESULTS 



The x^^ (zt) \/ (k^) extracted according to Eq. H22|l 
for various pxt and pTa. are shown in Fig. llfil and Fig. 1171 
In order to extract a (/c|,) values we have solved 



^hy(Pout>-(j|y)(4 + l)- ' (^t) = 

(49) 

for ylji^ where the (zt) and x\y—{pTa)/{pTt) are eval- 
uated according Eq. l|T7jl and Eq. (|^ respectively. 
These two quantities depend on (fc^) so we solved 
Eq. (|49|) iteratively by varying a (fc|,) value and in 
every step the (zt) and aih were recalculated. To do so 
the we need to know unsmeared final state parton spec- 
trum Tiq^pT) and the fragmentation function. For the 
latter one we used the LEP data (see Fig. 122(1 where the 
fragmentation functions of gluon and quark jets were 
measured in e+e" collision at y^=180 GeV. We have 
chosen 



Dl cx z-"(l-z)''(H-z)-T 



(50) 
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form used e.g. in 46] and extracted a, fi and 7 param- 
eters from the fit to distributions shown in Fig. (see 
Tab-EIIl. 



TABLE VIII: Extracted values of D[z) parameters accord- 
ing Eq. I|5()|l from the fit to the LEP data and power n of the 
unsmeared final state parton spectra T,q(pT) extracted from 
the fit to the single inclusive tt" invariant cross section (37l | 
for corresponding fragmentation function and fixed value of 
y{fcfy=2.5 GeV/c. 





gluon 


quark 


(gluon+quark)/2 


a 


0.16 


0.49 


0.32 




0.88 


0.57 


0.72 


7 


13.29 


8.00 


10.65 


n 


7.53 


7.28 


7.40 
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FIG. 26: (fey) values corresponding to Fig. 1161 as a solu- 
tion to Eq. (gnj for trigger tt" in 3 < prt < 4 GeV/c (sohd 
symbols) and 5 < prt < 10 GeV/c (open symbols) range. 
The solid and dashed lines bracket the systematic uncer- 
tainty due to the unknown ratio of quark and gluon jets, for 
the solid and open symbols, respectively. 

For a given set of parameters a, P and 7 the power 
of the unsmeared final state parton spectra was 
evaluated from the fit formula Eg . (I27II to the single 
inclusive tt'^ invariant cross section [33 ■ Here we used 
the simplified kx smearing 



Tt 



PTt 



-KiP' 



'Tt 



PT 



)exp ■ 



and for the fixed value of -y/ (k^) — \/2\/ C^tx) ~ ^-^ 
GeV/c the power n of Eg(pT) distribution was deter- 
mined. 

The measurement of the fragmentation functions at 
LEP was done separately for quark and gluon jets and 



FIG. 27: (fcy) values corresponding to Fig. 1171 as a so- 
lution to Eq. 1491 for associated particles in 1.4 < pTa < 
5 GeV/c region (solid symbols). The solid lines bracket the 
systematic error due to the unknown ratio of quark and 
gluon jets. The CCOR measurement at ^/s=Q2.4: GeV Q 
(empty triangles). 



the slopes of these two D{z) distributions are different. 
Quark jets produce a significantly harder spectrum than 
gluon jets (see Fig.l22|l. Since the relative abundance of 
quark and gluon jets at a/s=200 GeV is not known, for 
the final results we assumed that the numbers of quark 
and gluon jets are equal; the final D(z) uses the av- 
eraged parameter values between quark and gluon and 
the difference was used as a measure of the systematic 
uncertainty. 

Resulting a/ (fcy) values for 3 < prt < 4 GeV/c and 
5 < PTt < 10 GeV/c as a function of pTa are shown 
in Fig. [2E1 (compare to uncorrected values Fig. ITTjl . 
The solid and dashed lines bracket the systematic er- 
ror due to the unknown ratio of quark and gluon 
jets. These data points correspond to the uncorrected 
^h^ i^t) \/ (^t) values shown in Fig. El The \/Jk^ 
values for varying pTt corresponding to the data shown 
of Fig.Elare shown in Fig.|57| Also here the solid lines 
bracket the systematic error due to the unknown ra- 
ti o of q uark and gluon jets. It is evident that unfolded 
values reveal, within the error bars, no depen- 
dence neither on nor on pTt- The tabulated data 
are given in Table Hxl 

We compared the y/ data obtained in this anal- 
ysis to a/ {k"^) values found by the CCOR collaboration 
at Vs=62.4 GeV Q (empty triangles on Fig.EZJ- Al- 
though the trend with pTt seems to be similar the over- 
all magnitude at y/s=200 GeV is significantly higher. 

The (zt) and values from the iterative solution of 
Eq. (|49|) as a function of the Tr^trigger momenta p^t 
and associated momenta pTa are shown in Fig. 1281 and 
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FIG. 28; (^t) and as a function of prt for the 1.4 < 
PTa < 5.0 GeV/c associated region. 



FIG. 29: The (zt) and ih values (see Eg. 1201 as solution of 
Eq. (BHJ for 3 < pn < 4 GeV/c and 5 < pTt < 10 GeV/c 
as a function of pxa. ■ 



TABLE IX: Values of x'^ (zt) ^J {k^) and ^{fcf) for var- 
ious trigger particle prt and associated momenta in the 
1.4 < pTa < 5.0 GeV/c region shown in Fig. Eland Fig. 1771 



PTt 

GeV/e 



GeV/e 



GeV/c 



3.22 
3.89 
4.90 
5.91 
7.24 
8.34 



1.63 ± 
1.66 ± 
1.89 ± 
2.06 ± 
2.17 ± 
2.53 ± 



0.08 
0.08 
0.13 
0.19 
0.25 
0.62 



2.79 ± 
2.57 ± 
2.66 ± 
2.74 ± 
2.83 ± 
3.11 ± 



0.13 ± 
0.11 ± 
0.17 ± 
0.20 ± 
0.25 ± 
0.60 ± 



0.35 
0.33 
0.35 
0.34 
0.32 
0.33 



Fig. 123 There is an opposite trend; whereas the (zt) 
rises with prt it is falhng with pre,- It is an interesting 
consequence of two effects: competition between steeply 
falHng final state parton spectra and rising fragmenta- 
tion function with parton momentum. Secondly, the 
detection of trigger particle biases the kx vector in the 
direction of the trigger jet as discussed in section IVIBI 
The PTt averaged value of \J (k^) (Fig. \T7\i is com- 
pared to the average parton pair momentum, (p^) 
= {PT)pair, presented in 0| (see Fig. 1201 . The value 
of {pT)pair is determined as a sum of the two par- 
tons' (/ct)- In the present analysis the \/{k^ is de- 
termined and thus the value of {pT)pair is evaluated as 

iPT}pair = V2X {kr) = v/V2 X ^/W) ■ 

The present value of {pT)pair 



{P7 



' pair 



= 3.36 ± 0.09(stat) ± 0.43(sys) GeV/c 



appears to be in a good agreement with the lower energy 
dijet and dilepton measurements or the higher energy 
measurement in diphoton production . A UA2 mea- 
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FIG. 30: (color online) Compilation of mean pair pr mea- 
surements |2^ and comparisons to the {pT)pair measured in 
this analysis. 



surement of (pr) of Z° production at ^/s ^ 600 GeV 
gives 8.6 ± 1.5 GeV/c Uli^. 



VIII. SUMMARY 

We have made the first measurement of jet jr and 
kT for p -I- p collisions at ^/s = 200 GeV using the 
method of two-particle correlations. Analysis of the 
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TABLE X: The x"^ (zt) ^ {k-D yj {k-D values as a function of pTa for two different trigger tt transverse momentum bins 
shown in Fig. 1161 and Fig. 1261 AU units in rad and GeV/c. 



PTa 


o pi t ^ ^ 




PTa 


o \ pi t ^ J-U 


^{kl) 


1.7 


1.76 ± 0.12 


2.66 ± 0.19 


1.9 


2.69 ± 0.37 


3.09 ± 0.30 


2.2 


1.74 ± 0.13 


2.94 ± 0.22 


2.2 


2.54 ± 0.31 


3.19 ± 0.30 


2.7 


1.37 ± 0.13 


2.57 ± 0.23 


2.7 


2.13 ± 0.26 


3.04 ± 0.30 


3.2 


1.45 ± 0.12 


2.93 ± 0.23 


3.4 


1.89 ± 0.27 


3.04 ± 0.38 


3.9 


1.44 ± 0.11 


3.19 ± 0.23 


4.7 


1.41 ± 0.30 


2.64 ± 0.56 


5.0 


1.04 ± 0.10 


2.68 ± 0.25 









error. The final extracted values of kr are then de- 
Fl^m ^"^^ """"^ "^'^^ termined to be also independent of the transverse mo- 

ig- 1— Qj mentum of the trigger tt", in the range measured, with 

PTt ( GeV/c) (2t) values of .J {k"^) = 2.68 ± 0.07(stat) ± 0.34(sys) GeV/c. 

3!22 0.51 ± 4.10"'' ± 0.06 0.88 ± 0.01 

3.89 0.56 ± 2.10"^ ± 0.07 0.87 ± 0.01 

4.90 0.61 ± 1.10"^ ± 0.07 0.85 ± 0.01 

5.91 0.64 ± 1.10"'' ± 0.07 0.85 ± 0.02 
7.24 0.66 ± 1.10"^ ± 0.07 0.86 ± 0.02 

8.34 0.68 ± 5.10"=* ± 0.06 0.84 ± 0.05 Acknowledgments 



TABLE XII: The (zt) and ih values with pTa for two trigger 
tt" momenta bins as shown on Fig. 1291 

3 < PTt < 4 GeV/c 



PTa 




*h 


1.72 


0.54 ± 8.10"^ ± 0.06 


0.81 ± 0.01 


2.22 


0.52 ± 6.10"^ ± 0.06 


0.88 ± 0.01 


2.73 


0.51 ± 1.10"^ ± 0.07 


0.95 ± 0.01 


3.23 


0.49 ± 1.10"^ ± 0.06 


0.99 ± 0.01 


3.93 


0.47 ± 5.10"^ ± 0.06 


1.04 ± 0.01 


5.04 


0.41 ± 6.10"^ ± 0.06 


1.06 ± 0.01 




5 < PTt < 10 GeV/c 




PTa 






1.85 


0.66 ± 4.10"^ ± 0.06 


0.75 ± 0.04 


2.24 


0.64 ± 1.10~^ ± 0.06 


0.80 ± 0.03 


2.73 


0.61 ± 2.10"^ ± 0.07 


0.87 ± 0.02 


3.44 


0.57 ± 2.10"^ ± 0.07 


0.92 ± 0.02 


4.65 


0.52 ± 5.10"^ ± 0.08 


0.98 ± 0.01 



angular widths of the near-side peak in the correla- 
tion function has determined that the jet fragmenta- 
tion transverse momentum jt is constant with trigger 
particle pTt and the extracted value yj = 585 ± 
6(stat) ± 15(sys) MeV/c is comparable with previous 
lower i/i measurements. The width of the away-side 
peak is shown to be a measure of the convolution of 
jx with the jet momentum fraction z and the partonic 
transverse momentum /ct- (^t) is determined through 
a combined analysis of the measured tt*' inclusive and 
associated spectra using the jet fragmentation functions 
from e^e~ measurements. The average of (zt) from the 
gluon and quark fragmentation functions is used and 
the difference is taken as the measure of the systematic 
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APPENDIX A 

First and second moments of normally 
distributed quantities 



comparison of pxa. • (f^ ^ '^a + I'-'^A • ■ ■) with the true 
Pout magnitude (simple monte carlo) for various cta 



1.1r 



Let a; be a ID variable with normal (Gaussian) distri- 
bution and r — \J -V 'ip' is a 2D variable with x and y 
of normal distribution then the following relations can 










(r) 






(N> 






(kl> 




(r) 






-I 






2^2 = al 



Both and fcy are two dimensional vectors. We 
assume Gaussian distributed x and y components and 
thus the mean value {Rtx) and {kTy) is equal to zero. 
The non-zero moments of 2D Gaussian distribution are 
e.g. the root mean squares -\/ (j^), \/ {k'^) or the mean 
absolute values of the jt, projections into the per- 
pendicular plane to the jet axes (IjTyl) and (Ifcryl). 
Note that there are a trivial correspondences 

^)^^{kT) = V^{\kTy\) (Al) 



The correct way to analyze the azimuthal 
correlation function. 



V 0.9- 

3 0.8- Taylor expansion 

cvjj PjaSino^ 

^ 0.7- Pia^A 

- exp(-sinV2<p^^|))-cos((|)) 

0.& 

0.1 0.2 0.3 0.4 0.5 



0.6 0.7 



FIG. 31: (color online) The relative error on pout de- 
termination from the azimuthal correlation function based 
on the Taylor expansion of ^sinA<^^^ (dashed line), with 

an assumption of {Pout}~PTs. sincA (dotted line) and 

\/ (Pout)— PTaCTA (dotted-dashed line). The solid red line 

corresponds to \/ (Pout) from Eq. 1231 . 



Construction and fitting of the two-particle azimuthal 
correlation function is discussed in section HVl Tradi- 
tionally the correlation function is fitted by two Gaus- 
sian functions - one for intra-jet correlation (near peak) 
and one for the inter-jet correlations (away-side peak). 
From the extracted variances of the Gaussian functions 
the Jt and pout magnitudes are extracted. 

There is, however, a fundamental problem with this 
approach. The pout-vector defined in Eq. H17|l is equal 
to pTa sin A(j) event by event. However, we measure the 
width of the correlation peak and this corresponds to 
a/ (A(/)2) = a A- The relation \/ (Pout)~ PTaSinaA is 
not a good approximation for (ta>0.4 rad (see Fig. l^ . 
The assumption that the away-side correlation has a 
Gaussian shape is also good only for small values of cta 
(see Fig. ED). 

One way of relating VvPout) and a a was proposed e.g. 
by Peter Levai ^3 and used in several other analyzes. 



Since V {vlnt)=PT. 
relate Pout and a a is to expand 



^^ (siv? l^fjyj one possibility how to 



^sin^ A(? 



_ iA04 



45 



4 



+ gC'A ■ 



where we assumed a Gaussian distribution of Ac/). The 



values is shown in Fig. |^ It is obvious that there 
is only a little difference between VvPoiit) ^-P^a sintTA, 
V (Pout)=PTaO'A and the Taylor series. In the region 
where (ta>0.4 rad, all approximations seems to be 
equally bad. 

However, Pout, the only quantity with a truly Gaus- 
sian distribution (if we neglect the radiative corrections 
responsible for non-Gaussian tails in the Pout distribu- 
tion which are anyway not relevant for the kx analysis) 
can be directly extracted from the correlation function. 
With the assumption of Gaussian distribution in Pout, 
we can write the away-side A(/i-distribution (normalized 
to unity) as 



way 



dAcj) 



37r/2 
ir/2 



_ dN dpout 
dpout dAcj) 
Pts, cos A4) 



v/M^Erf(^^ 



) 



■ exp 



p^^ sin^ A0 
2(p2^,) 



This is the correct way of extracting a dimensional 
quantity from the azimuthal correlation function in the 
case of narrow associated bin. Similar line of arguments 
can be drawn also in the case of near peak. However, 
given the narrowness of the near angle peak, the simple 
Gaussian approximation is good enough. 



25 



[1] A. L. S. Angelis et al. (CERN-Columbia-Oxford- 
Rockefeller), Phys. Lett. B97, 163 (1980). 

P. Darriulat ct al., Nucl. Phys. B107, 429 (1976). [26 
M. Delia Ncgra et al. (CERN-Collogo do Francc- 
Heidelberg-Karlsruhe), Nucl. Phys. B127, 1 (1977). [27" 
K. Adcox et al. (PHENIX), Phys. Rev. Lett. 88, 022301 

(2002) , nucl-ex/0109003. [28 
S. S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, [29^ 
072303 (2003), nucl-ex/0306021. 

S. S. Adlcr ct al. (PHENIX), Phys. Rev. Lett. 91, [30 

182301 (2003), nucl-cx/0305013. 

C. Adlcr ct al. (STAR), Phys. Rev. Lett. 90, 032301 [31 

(2003) , nucl-ex/0206006. [32 
C. Adler et al. (STAR), Phys. Rev. Lett. 90, 082302 [33 

(2003), nucl-ex/02 10033. 

A. B. Migdal, Phys. Rev. 103, 1811 (1956). [34 
X.-N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 1480 
(1992). [35 
X.-N. Wang, Phys. Rev. C58, 2321 (1998), hep- 
ph/9804357. [36; 
C. A. Salgado and U. A. Wiedemann, Phys. Rev. Lett. 
93, 042301 (2004), hcp-ph/0310079. [37" 
J.-w. Qiu and I. Vitcv, Phys. Lett. B570, 161 (2003), 
nucl-th/0306039. [38 
X.-N. Wang, Nucl. Phys. A702, 238 (2002), hep- 
ph/0208094. [39; 
S. M. Herman, J. D. Hjorken, and J. B. Kogut, Phys. 
Rev. D4, 3388 (1971). [40 
J. F. Owens and J. D. Kimel, Phys. Rev. D18, 3313 
(1978). [41 
J. F. Owens, E. Reya, and M. Gluck, Phys. Rev. D18, 
1501 (1978). [42 
R. P. Feynman, R. D. Field, and G. C. Fox, Phys. Rev. 
D18, 3320 (1978). [43 
J. F. Owens, Rev. Mod. Phys. 59, 465 (1987). 
G. Bunce, N. Saito, J. Soffer, and W. Vogelsang, Ann. [44' 
Rev. Nucl. Part. Sci. 50, 525 (2000), hep-ph/0007218. [45 
R. Cutler and D. W. Sivers, Phys. Rev. D17, 196 
(1978). [46 
R. Cutler and D. W. Sivers, Phys. Rev. D16, 679 
(1977). [47; 

B. L. Combridge, J. Kripfganz, and J. Ranft, Phys. [48 
Lett. B70, 234 (1977). [49 
R. P. Feynman, R. D. Field, and G. C. Fox, Nucl. Phys. [50 
B128, 1 (1977). 

Y. L. Dokshitzer, V. A. Khoze, A. H. Mueller, and 



S. I. Troian, Basics of perturbative qcd, gif-sur-Yvette, 

France: Ed. Frontieres (1991) 274 p. (Basics of). 

L. Apanasovich et al., Phys. Rev. D59, 074007 (1999), 

hep-ph/9808467. 

A. Kulesza, G. Sterrrian, and W. Vogelsang, Nucl. Phys. 
A721, 591 (2003), hep-ph/0302121. 

F. W. Busser et al., Phys. Lett. B46, 471 (1973). 

S. M. Berman, J. D. Bjorken, and J. B. Kogut, Phys. 
Rev. D4, 3388 (1971). 

R. Blankenbecler, S. J. Brodsky, and J. F. Gunion, 
Phys. Lett. B42, 461 (1972). 

D. Antreasyan et al., Phys. Rev. D19, 764 (1979). 
P. Darriulat, Ann. Rev. Nucl. Part. Sci. 30, 159 (1980). 
K. Adachi et al. (TOPAZ), Phys. Lett. B451, 256 

(1999) , hep-ex/9901036. 

K. Adcox et al. (PHENIX), Nucl. Instrum. Meth. 
A499, 469 (2003). 

L. Aphecetche et al. (PHENIX), Nucl. Instrum. Meth. 
A499, 521 (2003). 

K. Adcox et al. (PHENIX), Nucl. Instrum. Meth. 
A499, 489 (2003). 

S. S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 

241803 (2003), hep-ex/0304038. 

L. Aphecetche et al. (PHENIX), Nucl. Instrum. Meth. 
A499, 521 (2003). 

M. Aizawa et al. (PHENIX), Nucl. Instrum. Meth. 
A499, 508 (2003). 

S. S. Adler et al. (PHENIX), Phys. Rev. C69, 034910 
(2004), nucl-cx/0308006. 

S. S. Adler et al. (PHENIX), Phys. Rev. Lett. 94, 
082301 (2005), nucl-ex/0409028. 

P. Levai, G. Fai, and G. Papp, Dijet correlations at isr 

and rhic energies, hep-ph/0502238. 

J. Adams et al. (STAR), Phys. Rev. Lett. 91, 072304 

(2003), nucl-cx/0306024. 

G. W. van Apeldoorn et al., Nucl. Phys. B91, 1 (1975). 
X.-N. Wang, Phys. Lett. B595, 165 (2004), nucl- 
th/0305010. 

P. Abreu et al. (DELPHI), Eur. Phys. J. C13, 573 

(2000) . 

G. Alexander et al. (OPAL), Z. Phys. C69, 543 (1996). 
F. Abe et al. (CDF), Phys. Rev. Lett. 70, 2232 (1993). 
R. Ansari et al. (UA2), Phys. Lett. B194, 158 (1987). 
R. Ansari et al. (UA2), Z. Phys. C41, 395 (1988). 



Ta 



Pt 



,pair 




T,pair 




